A Quark- Antiquark Condensate in Three-Dimensional QCD 
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Three-dimensional lattice QCD is studied by Monte Carlo simulations within the quenched ap- 
proximation. At zero temperature a quark- antiquark condensate is observed in the limit of vanishing 
quark masses. The condensate vanishes continuously at the finite-temperature deconfinement phase 
transition of the theory. A natural interpretation of this phenomenon in the full theory with dynam- 
ical quarks is in terms of the spontaneous flavor symmetry breaking U(Nf) — > U(Nt/2) x U(Nf/2). 
In addition, the spectrum of low-lying Dirac operator eigenvalues is computed and found to be 
consistent with a flat distribution at zero temperature, in agreement with analytical predictions. 
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Although super-renormalizable, non-Abelian gauge theories in (2+1) dimensions are notoriously difficult 
to analyze in perturbation theory due to the strong infrared divergences that appear already at leading 
orders. For this reason, a very appropriate framework for the study of such three-dimensional theories 
is that of lattice gauge theory In this letter we shall report on a series of Monte Carlo simulations 
aimed at clarifying the issue of flavor symmetry breaking and, at high temperature, flavor symmetry 
restoration in three-dimensional QCD. 

One striking difference between massless (2+l)-dimensional QCD and its (3+l)-dimensional counter- 
part is the absence of chiral symmetry. The lowest representation of the Clifford algebra in three 
dimensions has dimension two. Spinors are thus naturally taken to be two-spinors, and as 7-matrices 
one can choose the three Pauli matrices tTj. On the surface there seems to be no room for an additional 
analogue of 75, and hence no chiral symmetry. A fermionic mass term does, however, break parity V 
and time-reversal T symmetries 0, symmetries that are both respected by the massless Lagrangian. 
By introducing an even number of quark species Nf, and combining the spinors in pairs into four- 
spinors ip, one can rewrite the fermionic part of the Lagrangian entirely in terms of three 4-dimensional 
7-matrices 70,71,72 an d the four-spinors ip. A mass term of the kind rmpi/j can now be introduced 
H], |3). Such a mass term is V and T invariant. It does, however, break both of the "chiral" symmetries 
associated with 74 and 75 rotations. If we restrict our attention to an even number of massless quark 
flavors we thus have a formulation which is very similar to that of four-dimensional massless QCD. 
There are chiral symmetries, and the strong, in the absence of matter fields almost certainly confining, 
gauge interactions. For a small even number of flavors one could expect these chiral symmetries to 
break spontaneously. 

It has been suggested that a plausible flavor symmetry breaking in the theory with an even number of 
flavors Nf should be that of U(N f ) -> U(N f /2) x U(N f /2) [|, f§. The formation of a chiral condensate 
(ijjip) would be a signal for such spontaneous symmetry breaking. If correct, it would open up the 
possibility of comparing the detailed universality predictions [|], ||] for the Dirac operator spectrum 
around eigenvalues near A = with lattice Monte Carlo data, as has recently been done in the case of 
four-dimensional lattice QCD H. The first step in such a program is to establish the existence of the 
chiral condensate {tfiiff). 

The transcription of the formulation in terms of four-spinors and four-dimensional 7-matrices to a 
lattice theory is quite naturally done in terms of staggered fermions p[|. Interpreting three-dimensional 
staggered fermions in this manner has already been explored in several Monte Carlo analyses, see e.g. 
ref. ||. The standard three-dimensional action for free staggered fermions, keeping the lattice spacing 
a explicitly, 

S F + S M = ~\a d - 1 J2(-l) n+r2+ - +r ^ [x(r)x(r + /x) + x(r + /x)x(r)] + ima d £ x(r) X (r) (1) 

r,fi r 

can be transformed into a form with 2-component Dirac fields. In the limit a— >0 it will go to the 
continuum action for free Dirac fermions with two flavors and (if the kinetic terms is defined with 
respect to the same 7 matrices; see below) masses m and —m. In three dimensions there are two 
equivalence classes of irreducible representations of the Clifford algebra, {7^} and {— 7^} = {0^}. 

Following ref. 0] the new action can be defined on a lattice with twice the original lattice spacing: 
S F + S M = (2a d ) J2 {"(7/, ® I)dpU + 303„ ® I)dpd + a[u(I ® 7 J)Sjd + 3(1 ® P^u] 

+ im[u(I® I)u + d(I ® I)d\ \ (2) 



2 



Here u and d are 2-component spinors and in the quark bilinears the first 2x2 matrix acts on spinor 
indices while the second acts on flavor indices. One can identify a U(l) x U(l) symmetry in eq. (Q) as 

(u,d) -> (u,d)e- m (3) 



(u,d)^(u,d)[.. Q . (4) 







cos 02 i sin 62 
i sin 62 cos #2 

The symmetry in eq. (Q) is a remnant of the three-dimensional "chiral symmetry" on the lattice, and 
it is broken by the inclusion of a mass term. The symmetry of eq. (||) is nothing but fermion number 
conservation. 

We consider gauge group SU(S), and add to the fermionic part of the action the conventional Wilson 
gauge action (/? = 6/g 2 ): 

Sgauge = /?E( 1 -^ ReTrC/ p) . ( 5 ) 
p 6 

where U p is the product of links along a fundamental plaquette. The canonical ensemble of gauge 
field configurations was generated using a combination of microcanonical over-relaxation (MOR) || 
and quasi-heat-bath (QHB) algorithms An individual MOR or QHB step consisted of updating 
consecutively the three 2x2 submatrices of the SU(3) link variable [|TTJ. In our implementation five 
MOR steps were followed by two QHB ones. 

We now present the results of our Monte Carlo simulations. Instead of first considering symmetric 
(2+l)-dimensional lattice volumes, we have immediately turned to lattices with the temporal direction 
much smaller than the two spatial directions. This has allowed us to study the finite-temperature 
behavior of the theory, while at the same time determining, in the low-temperature region, whether 
at all a quark-antiquark condensate is formed at zero temperature. As mentioned in the introduction, 
all simulations have been performed in the quenched approximation. This should have no bearing 
on the qualitative features as compared to the full theory with dynamical quarks, at least as long as 
the number of flavors is relatively small. We have monitored all relevant gauge-invariant correlation 
functions in both the matter and gauge sector, and in particular of course the crucial order parameters: 
the Polyakov line (W) and the chiral condensate {'ipip). Since the gauge dynamics is unaffected by the 
fermions in this approximation, the theory should undergo a 2nd-order finite-temperature deconfining 



phase transition at an intermediate 0- value [12]. This particular detail may of course be modified by 
the presence of (light) dynamical quarks (just as its four-dimensional counterpart), but the qualitative 
feature of a switch from confining to deconfined characteristics should persist. 

The finite-temperature phase transition in the gauge sector is associated to global Z(3) symmetry 
breaking. The Polyakov line aligns, in the deconfined phase, along the three complex roots of = 1 
on the unit circle, i.e along z £ {1, exp[2Mr/3], exp[4i7r/3]}. In the theory with dynamical quarks, these 
act like a magnetic field that force the spontaneous symmetry breaking on the Polyakov line to occur 
on the physical real axis (i.e. along z = l) fj~3| , but in this quenched situation there is equal probability 
for the Polyakov line to fall into any of the three sectors. This is normally not a question of great 
concern, since in finite-volume simulations one in any case needs a modified working definition of the 
Polyakov line expectation value, which customarily is based on taking its absolute value. But as has 
recently been emphasized [14], the existence of complex phases in the Polyakov line condensate can 



have strong effects on the chiral condensate in quenched simulations. f\ We have seen this very clearly 

1 The analogous phenomenon for SU(2) (where it just amounts to a sign, and in fact can effectively flip the fermion 
boundary conditions from antiperiodic to periodic), was already pointed out in ref. []15|. 
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Figure 1: The (3 dependence of the Polyakov line (W), measured on a 40 2 x 4 lattice. 

in our simulations. In Figure [j] we display the (absolute value) of the Polyakov loop expectation value 
on a 40 2 x 4 lattice, as a function of the gauge coupling (3. It shows the expected behavior of a smooth 
phase transition at an intermediate value of (3, which we read off to be around (3 ~ 15. This is in 
good agreement with earlier Monte Carlo results |l6| and mean field theory |j"7j| . It is also in rough 
agreement with three-dimensional scaling, and earlier Monte Carlo results for a temporal extent of 



N T = 2 which gave a critical /3-value of /3 C ~ 8.1 [12|. If we re-instate the lattice spacing a, the definition 
of (3 actually becomes f3 = 6/(g 2 a), that is, the coupling g is dimensionful in three dimensions. With 
temperature given by T = l/(N T a), perfect scaling in this theory would imply that for the physical 
critical temperature T c to stay fixed, the critical coupling f3 c should be 2 x 8.1 = 16.2 on our lattice with 
N T = 4. Thus, although continuum scaling is not perfectly obeyed here, the deviation is nevertheless 
only around 5%. 

While the behavior of the Polyakov line thus is completely as expected, the condensate {ipip) displayed 
at first a quite surprising behavior, which we eventually could attribute to the complex phases of the 
Polyakov line. Before presenting details on this particular issue, we first display, in Figure [2|, a typical 
plot of our final results (here for a quark mass of am q = 0.1). This figure, reminiscent of similar plots for 
SU(2) in 4 dimensions p5| , suggests immediately that, in the limit of massless quarks, (1) a condensate 
is formed in the zero-temperature theory (corresponding to the phase with (3 < f3 c ), but not at high 
temperature - we shall confirm these observations with explicit extrapolations to the massless limit 
below - and (2) at the deconfinement transition this condensate approaches zero smoothly. Neither 
of these observations come as a surprise. First of all, the formation of a -(/^-condensate at strong 
coupling is almost built in by the present lattice formulation. The strong-coupling, zero-temperature 



analysis of ref. [18] straightforwardly carries over to the present (2+l)-dimensional case, where it again 
predicts the formation of a chiral condensate.0 Variational calculations in the Hamiltonian formalism 
|j~9[ ] indicate the existence of the condensate also at finite couplings. Similarly, the extension to finite 
temperature |2(| immediately predicts, at strong coupling, a symmetry-restoring phase transition at 
finite temperature. 

2 The prediction is not as strong as in (3+1) dimensions, however, since one must view it in terms of a l/d-expansion 
that quantitatively could fail for d = 3. 
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Figure 2: The (3 dependence of the quark-antiquark condensate (ipifi), for am q = 0.1, measured on a 
40 2 x 4 lattice. 



Of course, details will depend on the quenched approximation. For instance, it is almost impossible to 
imagine that in the presence of a continuous deconfining phase transition the behavior of {ifii/j) could 
be anything but smooth in our case. To show that this is not what is found if one blindly measures 
(Tpi/j) in this quenched simulation, we show in Figure ^ (top) a typical run history of that quantity. 
For illustrative purposes we have here chosen /3 = 14.6 associated with Z(3) breaking. Surprisingly, the 
chiral condensate displays what appears to be a clear-cut two-state signal, indicating the proximity of 
a discontinuous phase transition. The reason for this bizarre behavior becomes evident if we compare 
(see Figure ^ (middle)) with the corresponding time history of the (negative of the) real part of the 
Polyakov line, (- ReVF). The apparent "two-state signal" appears to be reproduced in this variable, 
and in fact the data in the upper and middle part of Figure [3| are almost perfectly correlated. We recall 
that in the case of the Polyakov line, this behavior simply reflects the flipping around between the three 
different Z{3) phases, and in no way implies a two-state signal in the properly defined Polyakov line. 
Furthermore, the anti-correlated behavior between (ReW) and (ijjtp) is in complete agreement with the 
observations of ref . |l4j] : When the Polyakov line lies along one of its two complex phases (and its real 
part hence is relatively small) , the chiral condensate has an unphysical larger value than in the physical 
situation corresponding to a Polyakov line aligned along the real axis. 

The problem described above can be solved in different ways. We have opted for the following simple 
method. The phase eft of (W) is measured. Then, if eft > 7r/3, all the timelike link variables in the t = 
time slice are multiplied by exp(— 2ni/3). If 4> < — vr /3, these variables are multiplied by exp(27ri/3). If 
we find — 7r/3 < (p < vr/3, the configuration is left unchanged. Thus, for the transformed configuration, 
the phase of (W) lies between —tt/3 and 7r/3. The transformation as described is an exact symmetry of 
the quenched theory. Nevertheless, we prefer to discard the transformed configuration after performing 
measurements, and to continue the Monte Carlo process with the original configuration. In the bottom 
part of Figure ^ we show how the time history of (V?^) looks with this prescription. We have here used 
precisely the same starting configuration, and yet the result is radically different. There is now no trace 
of any metastability, and tpip is nicely fluctuating around one particular, stable, value. We have applied 
this prescription wherever we detected a two-state signal in the (tpip) time history. This resulted in a 
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Figure 3: The Monte Carlo time history, averaged over 5 consecutive iterations, of the quark-antiquark 
condensate (ipip) for am q = 0.1 measured on the original configurations (top), of the negative of the real 
part of (W) (middle) and of (ipip) measured on the transformed configurations (bottom) at (3 = 14.6 
on a 40 2 x 4 lattice. 
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Figure 4: The quark-antiquark condensate {ipip) as a function of the quark mass am q at (3 = 12.0, 
corresponding to the low-temperature phase, for various volumes. Also shown is the value for the 
condensate at am q = in the infinite volume limit (square). 



smooth dependence of (ipip) on ft, as shown in Figure |||] 

While the data already shown indicate the presence of a non- vanishing quark-antiquark condensate at 
low temperature, we have performed a much more detailed analysis to really confirm this result. The 
problem of course is that for any finite volume V there is no spontaneous symmetry breaking in the 
am q — > limit. Fortunately the behavior with increasing volume is, however, quite different in a phase 
of broken symmetry compared with a phase of restored symmetry. We can use this different behavior 
to test the hypothesis of spontaneous symmetry breaking in the confined phase of the theory. In Figure 
|2| we display data {ipip) for a variety of lattice volumes at = 12.0 (corresponding to the confined 
phase). While the condensate eventually goes towards zero as am q — * for any of the lattices volumes, 
there is a clear trend towards increasing values, at fixed am q , as the volume is increased. On the same 
figure we have also indicated the value of the condensate at am q = and extrapolated to V = oo 
as obtained from an entirely different analysis (see below). Considering the errors, this extrapolated 
value is completely consistent with the direct measurements at increasingly larger volumes, as shown 
in Figure It finally remains to be tested if the behavior on the high-temperature side is consistent 
with a vanishing condensate in the limit of zero quark masses. In Figure |5| we show the convergence 
of (V>V ; ) at /? = 17.0. These data are clearly very different from those in the confined phase, with no 
discernible volume dependence. All points lie on a straight line that directly extrapolates to a vanishing 
condensate at am q = 0. 

3 The used prescription is, however, not entirely free of systematic error: to every transition between the physical and 
the unphysical values of ^ip in Figure ^ corresponds a spike in the refined value of tpip. As a result of these spikes, we 
expect the average {4>ip) to be slightly biased upwards. This bias can be systematically reduced by increasing the lattice 
volume, thereby making the transitions between the phases of (W) increasingly rare. 
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Figure 5: Same as Figure |] but in the high-temperature phase, at /3 = 17.0. 



We have shown that, at least in the quenched approximation and in our present latticized version, a 
chiral condensate is formed in (2+l)-dimensional QCD. This condensate disappears according to expec- 
tations at the finite temperature deconfinement phase transition at T c . In the low-temperature phase 
it is thus possible to test the detailed predictions concerning the Dirac spectrum at small eigenvalues 
0, En. In contrast to the situation in four-dimensional QCD, the prediction is here that the quenched 
microscopic spectral density should be completely flat.0 To test this, we have performed a series of 
measurements of the lowest-lying Dirac operator eigenvalues, using the Ritz functional algorithm of ref. 
[21]. The density of low-lying eigenvalues at (3 = 12.0 is shown in Figure for lattice sizes 8 3 to 14 3 . 
Except for the largest lattice we computed the 10 lowest eigenvalues, on 14 3 the 16 lowest ones. There 
is a clear trend towards a flat spectral distribution with increasing lattice volume. The dip very near 
A ~ is due to finite-size effects, and indeed this dip becomes more and more narrow as the volume is 
increased. While a small peak appears close to A = 0, its significance is doubtful since it seems to be 
smallest on the largest lattice volume. In addition, it is not well correlated with the distribution of the 
smallest eigenvalue, as one would have expected from any genuine oscillatory behavior. From the flat 
plateaux we can extract an average value of the spectral density near the origin, i.e. p(0). We find a 
slight volume dependence in the extracted values of p(0). Extrapolating as 1/V to the infinite- volume 
limit we obtain ^(0) = 0.0315(19). Through a 3-dimensional analogue of the Banks-Casher formula this 
provides us with an independent measurement of the condensate: (ipip) = 7rp(0). Inserting the above 
value for p(0), this gives ("ijiip) = 0.099(6), which is the value shown in Figure It is consistent with our 
direct measurements. We have thus provided yet more independent confirmation of the formation of a 
condensate in the confined phase, while simultaneously testing the predictions about the microscopic 
spectral density of the Dirac operator [0, ||. 



4 This assumes that the quenched model with staggered quarks can be understood as the Nf 
with an even number of quarks. 
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Figure 6: The density of low lying eigenvalues for various volumes at (3 = 12.0. 
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